What Is Claimed Is: 



1 LA method for using a computer system to solve an unconstrained 

2 interval global optimization problem specified by a function / wherein/is a scalar 

3 function of a vector x = (x Ir x 2 , x 3> ... x n ), the method comprising: 

4 receiving a representation of the function / at the computer system; 

5 storing the representation in a memory within the computer system; and 

6 performing an interval global optimization process to compute guaranteed 

7 bounds on a globally minimum value of the function/x); 

8 wherein performing the interval global optimization process involves, 

9 applying term consistency over a subbox X, and 

1 0 excluding any portion of the subbox X that violates term 

1 1 consistency. 

1 2. The method of claim 1 , wherein applying term consistency 

2 involves: 

3 symbolically manipulating an equation within the computer system to 

4 solve for a first term, g(x'), thereby producing a modified equation g(x') = h(x), 

5 wherein the first term g(x') can be analytically inverted to produce an inverse 

6 function g'\y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X') = /i(X); 

9 solving for X ' = g ~\h(X)); and 

1 0 intersecting X ' with the subbox X to produce a new subbox X + ; 

1 1 wherein the new subbox X + contains all solutions of the equation within 

1 2 the subbox X, and wherein the size of the new subbox X + is less than or equal to 

1 3 the size of the subbox X. 
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1 3 . The method of claim 1 , wherein performing the interval global 

2 optimization process involves: 

3 keeping track of a least upper bound fjbar of the function/x); 

4 removing from consideration any subbox for which/x) > fjbar; 

5 applying term consistency to the inequality fix) < fjbar over the subbox X; 

6 and 

7 excluding any portion of the subbox X that violates the inequality. 

1 4. The method of claim 1 , wherein performing the interval global 

2 optimization process involves: 

3 determining a gradient g(x) of the function/x), wherein g(x) includes 

4 components g,(x) (z=l, ...,«); 

5 removing from consideration any subbox for which g(x) is bounded away 

6 from zero, thereby indicating that the subbox does not include a global minimum 

7 of/(x); and 

8 applying term consistency to each component g,(x)=0 (i=l, ...,/?) of g(x)=0 

9 over the subbox X; and 

10 excluding any portion of the subbox X that violates a component. 

1 5. The method of claim 1 , wherein performing the interval global 

2 optimization process involves: 

3 determining diagonal elements H n (x) (/=1, ...,«) of the Hessian of the 

4 function^*); 

5 removing from consideration any subbox for which a diagonal element of 

6 the Hessian is always negative, which indicates that the /is not convex and 

7 consequently does not contain a global minimum within the subbox; 
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8 applying term consistency to each inequality #„(x) > 0 (i=l, ...,«) over the 

9 subbox X; and 

1 0 excluding any portion of the subbox X that violates an inequality. 

1 6. The method of claim 1 , wherein performing the interval global 

2 optimization process involves performing the Newton method, wherein 

3 performing the Newton method involves : 

4 computing the Jacobian J(x,X) of the function / evaluated as a function of 

5 a point x over the subbox X; 

6 computing an approximate inverse B of the center of J(x ? X); and 

7 using the approximate inverse B to analytically determine the system 

8 Bg(x), wherein g(x) is the gradient of the function/x), and wherein g(x) includes 

9 components g t (x) (/=i, ...,«); and 

10 applying term consistency to each component (Bg (x)) ; = 0 (7=1, ...,ri) for 

1 1 each variable x t (/=1, ...,«) over the subbox X; and 

12 excluding any portion of the subbox X that violates a component. 

1 7. The method of claim 1 , further comprising terminating attempts to 

2 further reduce the subbox X when: 

3 the width of X is less than a first threshold value; and 

4 the magnitude of/(X) is less than a second threshold value. 

1 8. A computer-readable storage medium storing instructions that 



2 when executed by a computer cause the computer to perform a method for using a 

3 computer system to solve an unconstrained interval global optimization problem 

4 specified by a function/ wherein/is a scalar function of a vector x = {x }} x 2 , x 3 , ... 

5 x„) 9 the method comprising: 
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6 receiving a representation of the function/ at the computer system; 

7 storing the representation in a memory within the computer system; and 

8 performing an interval global optimization process to compute guaranteed 

9 bounds on a globally minimum value of the function fix); 

10 wherein performing the interval global optimization process involves, 

1 1 applying term consistency over a subbox X, and 

1 2 excluding any portion of the subbox X that violates term 

13 consistency. 

1 9. The computer-readable storage medium of claim 8 ? wherein 

2 applying term consistency involves: 

3 symbolically manipulating an equation within the computer system to 

4 solve for a first term, g(x') ? thereby producing a modified equation g(x ? ) = h(x\ 

5 wherein the first term g(x') can be analytically inverted to produce an inverse 

6 function g~ ; (y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X') = /?(X); 

9 solving for X ' = g~\h(X)); and 

1 0 intersecting X ' with the subbox X to produce a new subbox X + ; 

1 1 wherein the new subbox X + contains all solutions of the equation within 

12 the subbox X, and wherein the size of the new subbox X + is less than or equal to 

13 the size of the subbox X. 

1 10. The computer-readable storage medium of claim 8, wherein 

2 performing the interval global optimization process involves: 

3 keeping track of a least upper bound f_bar of the function/(x); 

4 removing from consideration any subbox for whichy(x) > fj>ar\ 
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and 



applying term consistency to the inequality fx) < fjbar over the subbox X; 
excluding any portion of the subbox X that violates the inequality. 



1 11. The computer-readable storage medium of claim 8, wherein 

2 performing the interval global optimization process involves: 

3 determining a gradient g(x) of the function/x), wherein g(x) includes 

4 components g,(x) (r= 1, . . . , n); 

5 removing from consideration any subbox for which g(x) is bounded away 

6 from zero, thereby indicating that the subbox does not include a global minimum 

7 of fx); an <3 

8 applying term consistency to each component g,(x)=0 (i=l, ...,«) of g(x)=0 

9 over the subbox X; and 

1 0 excluding any portion of the subbox X that violates a component. 

1 12. The computer-readable storage medium of claim 8, wherein 

2 performing the interval global optimization process involves: 

3 determining diagonal elements H n (x) (z = 1, . . . ,ri) of the Hessian of the 

4 function fix); 

5 removing from consideration any subbox for which a diagonal element of 

6 the Hessian is always negative, which indicates that the /is not convex and 

7 consequently does not contain a global minimum within the subbox; 

8 applying term consistency to each inequality H u (x) > 0 (7=1, ... s ri) over the 

9 subbox X; and 

1 0 excluding any portion of the subbox X that violates an inequality. 
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1 13. The computer-readable storage medium of claim 8 5 wherein 

2 performing the interval global optimization process involves performing the 

3 Newton method, wherein performing the Newton method involves: 

4 computing the Jacobian J(x ? X) of the function/ evaluated as a function of 

5 a point x over the subbox X; 

6 computing an approximate inverse B of the center of J(x ? X); and 

7 using the approximate inverse B to analytically determine the system 

8 Bg(x), wherein g(x) is the gradient of the function/x), and wherein g(x) includes 

9 components g,(x) ...,«); and 

1 0 applying term consistency to each component (Bg (x)X = 0 (i= 1, . . . , n) for 

1 1 each variable x, (i= 1, . . n) over the subbox X; and 

12 excluding any portion of the subbox X that violates a component. 

1 14. The computer-readable storage medium of claim 8, wherein the 

2 method further comprises terminating attempts to further reduce the subbox X 

3 when; 

4 the width of X is less than a first threshold value; and 

5 the magnitude of XX) is less than a second threshold value. 

1 1 5. An apparatus that solves an unconstrained interval global 

2 optimization problem specified by a function f wherein /is a scalar function of a 

3 vector x = (jc;, x% xs, ... x n \ the apparatus comprising: 

4 a receiving mechanism that is configured to receive a representation of the 

5 function/; 

6 a memory for storing the representation; and 
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7 an interval global optimization mechanism that is configured to perform 

8 an interval global optimization process to compute guaranteed bounds on a 

9 globally minimum value of the function /(x); 

10 a term consistency mechanism within the interval global optimization 

1 1 mechanism that is configured to, 

12 apply term consistency over a subbox X, and to 

1 3 exclude any portion of the subbox X that violates term 

14 consistency. 

1 16. The apparatus of claim 15 ? wherein the term consistency 

2 mechanism includes: 

3 a symbolic manipulation mechanism that is configured to symbolically 

4 manipulate an equation within the computer system to solve for a first term, g(x') ? 

5 thereby producing a modified equation g(x') = h(x), wherein the first term g(x') 

6 can be analytically inverted to produce an inverse function g'\y); 

7 a solving mechanism that is configured to, 

8 substitute the subbox X into the modified equation to 

9 produce the equation g(X 9 ) = h(X) y and to 

1 0 solve for X' = g } ^(h(X)); and 

1 1 an intersecting mechanism that is configured to intersect X' with the 

1 2 subbox X to produce a new subbox X + ; 

1 3 wherein the new subbox X + contains all solutions of the equation within 

14 the subbox X, and wherein the size of the new subbox X + is less than or equal to 

15 the size of the subbox X. 

1 17. The apparatus of claim 1 5, 

2 wherein the interval global optimization mechanism is configured to, 
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1 keep track of a least upper bound fjbar of the function/[x), 

2 and to 

3 remove from consideration any subbox for which 

4 /x) >f_bar; 

5 wherein the term consistency mechanism is configured to, 

6 apply term consistency to the inequality fix) < fjbar over 

7 the subbox X, and to 

8 exclude any portion of the subbox X that violates the 

9 inequality. 

1 18. The apparatus of claim 1 5 , 

2 wherein the interval global optimization mechanism is configured to, 

3 determine a gradient g(x) of the function/x), wherein g(x) 

4 includes components g,(x) (z=l, .„,/?), and to 

5 remove from consideration any subbox for which g(x) is 

6 bounded away from zero, thereby indicating that the subbox does 

7 not include a global minimum of/(x); and 

8 wherein the term consistency mechanism is configured to, 

9 apply term consistency to each component gi(x)=0 

10 (zW, of g(x)=0 over the subbox X, and to 

1 1 exclude any portion of the subbox X that violates a 

12 component. 

1 19. The apparatus of claim 1 5 , 

2 wherein the interval global optimization mechanism is configured to, 

3 determine diagonal elements H n (x) (7=1, ...,n) of the 

4 Hessian of the function/(x), and to 
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5 remove from consideration any subbox for which a 

6 diagonal element of the Hessian is always negative, which 

7 indicates that the /is not convex and consequently does not contain 

8 a global minimum within the subbox; and 

9 wherein the term consistency mechanism is configured to, 

10 apply term consistency to each inequality H u (x) > 0 

1 1 (z=l, ...,«) over the subbox X, and to 

12 exclude any portion of the subbox X that violates an 

13 inequality. 

1 20. The apparatus of claim 15, further comprising a Newton 

2 mechanism within the interval global optimization mechanism that is configured 

3 to: 

4 compute the Jacobian J(x ? X) of the function / evaluated as a function of a 

5 point x over the subbox X; 

6 compute an approximate inverse B of the center of J(x,X); and to 

7 using the approximate inverse B to analytically determine the system 

8 Bg(x) ? wherein g(x) is the gradient of the function/x), and wherein g(x) includes 

9 components g z ( x ) (z=l, ... } n); 

10 wherein the term consistency mechanism is configured to, 

1 1 apply term consistency to each component (Bg (x)) z = 0 

12 (zW, ...,«) for each variable x t ...,«) over the subbox X, and to 

1 3 exclude any portion of the subbox X that violates term 

14 consistency. 
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1 21. The apparatus of claim 15, further comprising a termination 

2 mechanism that is configured to terminate attempts to further reduce the subbox X 

3 when: 

4 the width of X is less than a first threshold value; and 

5 the magnitude of/(X) is less than a second threshold value. 
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